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We present a general theory for the full counting statistics of multiple Andreev reflections in 
incoherent superconducting-normal-supcrconducting contacts. The theory, based on a stochastic 
path integral approach, is applied to a superconductor-double barrier system. It is found that all 
cumulants of the current show a pronounced subharmonic gap structure at voltages V = 2A /en. For 
low voltages F < A/e, the counting statistics results from diffusion of multiple charges in energy 
space, giving the pth cumulant (Q p ) cx V 2_p , diverging for p > 3. We show that this low-voltage 
result holds for a large class of incoherent superconducting-normal-superconducting contacts. 
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Electrical transport in superconducting contacts is a 
subject of persistent interest. The basic mechanism of 
transport in voltage biased contacts, Multiple Andreev 
Reflections (MAR), was described by Klapwijk, Blonder 
and Tinkham []J twenty years ago. Since then, a large 
number of works have investigated the current-voltage 
characteristics in superconducting contacts. As a promi- 
nent example, experiments on superconducting atomic 
point contacts have been explained remarkably well by 
coherent MAR-theories Q. 

Recently, interest has turned to the properties of cur- 
rent noise. Shot noise was measured in several types 
of superconducting contacts 0, 0, 01 and was theoret- 
ically studied in both the coherent j^, and incoher- 
ent limit 0, 0, 0|. Very recently, the full distribu- 
tion of current fluctuations, the Full Counting Statistics 
(FCS) [m of MAR in coherent superconducting con- 
tacts was calculated 0, The FCS revealed ex- 
plicitly that charge in coherent supcrconducting-normal- 
supcrconducting (SNS) systems is transfered in multiples 
of the elementary charge. In diffusive SNS-systems, this 
leads to a characteristic power-law divergence of the shot 
noise at low voltages, Si ~ V^ 1 / 2 0. However in many 
experiments quantum coherence is suppressed |fj . Inter- 
estingly, calculations for incoherent diffusive structures 
Q EI [ID show that the low-voltage shot noise is well- 
behaved and can be explained in terms of an effective 
electron temperature of the order of the superconducting 
gap A. This naturally raises the question about the role 
of multiple charge transfer in incoherent SNS-systems. 

We address this question by developing a theory for the 
FCS of incoherent MAR in semiclassical SNS-systems, 
based on a stochastic path integral 0, ^| approach. 
As an illustrative example, we investigate the supercon- 
ducting double-barrier model of Octavio et al. (OTBK) 
|l6j| . focusing on the voltage dependence of the FCS. For 
a generic incoherent SNS-system, we show that for low 
voltage, eV <C A, the generating function is 



with t the measurement time, \ the counting field and G 
the linear, low-voltage conductance. This result, an effect 
of diffusive motion of the quasiparticles in energy space, is 
just the FCS for a metallic diffusive wire 0,03 with an 
effective voltage dependent charge e — > e(2A/eV). Un- 
like coherent diffusive systems [tJ, the cumulants (Q p ) = 
d p S/dx p \ x =v ~ V 2 ~p diverge only for p > 3. The suc- 
cessful measurement of the third cumulant (Q 3 ) in nor- 
mal conducting tunnel junctions shows that an ex- 
perimental test of our prediction is feasible. 
The model - The studied geometry is shown in Fig. 
H Two superconducting electrodes (S) with gaps A 
are connected via a ballistic normal conducting re- 
gion (N) with M transport modes. The NS-interfaces, 
with equal mode-independent normal state transparen- 
cies r, are characterized by energy dependent probabili- 
ties |l9j for normal/ Andreev reflection and transmission 
Rn/a{E) 1 T n / A {E). A voltage V is applied between the 
two superconductors. For eV < 2A, transport of quasi- 
particles across the gap takes place via MAR along the 
energy ladder shown in Fig. ^ 

In contrast to previous works on FCS of MAR fl^ll3j| . 
we will consider the incoherent limit, i.e. with suppressed 
proximity and dc/ac-Josephson effect. This experimen- 
tally interesting limit || is relevant for strong dephasing 




is 


E' 


1 




T 
V 








SNS,2n+l 


1 



RXE) 




s(x) = 



2A 



-asinhVexp(x2A/eU) - 1, (1) 



FIG. 1: Left panel: A schematic of the geometry and the 
MAR- ladder in energy space (see text for details). Right 
panel: Energy dependent probability Ra(E) of a single An- 
dreev reflection for different interface transparencies F. 
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in the normal region (for incoherent transport in diffu- 
sive NS-systems, see Refs. Hi})- We assume low temper- 
ature, kT -C A, and negligible relaxation due to electron- 
electron or electron-phonon scattering. Furthermore, we 
focus on semiclassical systems with M 3> 1. 

In this situation, following OTBK [lfj| . the state of 
the normal region is described by energy, space and di- 
rection dependent occupation functions /„, /„ for left- 
and right-going particles (n even) and holes (n odd) [see 
Fig. |IJ. The average current through the junction is ob- 
tained from a set of coupled linear equations, the OBTK- 
equations, for the occupation functions on different rungs 
of the energy ladder. In the present work, we show how 
to extend the OBTK-equations to determine the current 
noise and the FCS in the system. We however emphasize 
that incoherent diffusive [alilEl or chaotic contacts [2l| 
can be treated in a very similar fashion. Below we put 
e = ft = 1. 

Introduced by Levitov and Lesovik ^lj, the FCS 
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is the distribution giving the probability that a certain 
charge Q flows through an electrical conductor during 
the time interval [0, r]. It describes entirely the low- 
frequency current fluctuations, if r is chosen longer than 
any intrinsic time scale. The Fourier transform of the 
probability distribution P{Q) yields the generating func- 
tion S(x) = X p (Q p ) / P^- °f irreducible cumulants and 
gives the mean current I = (Q)/r, the current noise 
Si = (Q 2 )/t, the third cumulant C 3 = (Q a )/r describing 
the asymmetry of the distribution, etc. 

Formal Solution - The FCS of a single NS-interface is 
a multinomial process and can be fully described j2^| by 
the reflection and transmission probabilities R? 
Inside the gap, \E\ < A, the generating function is 



^ \n\fafa + f p (l - f h )(R A e 2x + Rn) 



+fa(l - f P )(R A e-^ + R N ) + (1 - / P )(l - fa)} (3) 

where /„, fa are the particle and hole occupation func- 
tions in the normal region. A similar expression accounts 
for reflection and transmission outside the gap, \E\ > A. 
The SNS-junction can be considered as a series of two 
NS-intcrfaces sharing the same occupation functions. In 
the incoherent semiclassical limit, the FCS of such a com- 
pound system can be calculated with the stochastic path 
integral formalism [Til fl5| . The key ingredient in this 
formalism is a separation of time scales for fast collision 
events (the Andreev reflections) and for the slow evolu- 
tion of the occupation functions. This separation allows 
one to integrate out the fast fluctuations and to express 
the generating function S(x) by an action in terms of 



slow collective variables as 
r 2V 



S\x,f,X) = I dE S NS ,n(E n ), (4) 
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where the integral is over the different ladders, contribut- 
ing independently to the action. The actions Sns,™ eval- 
uated at E n = E + nV (see Fig. ^) arc generalizations of 
Eq. which distinguish left- and right-going currents. 
E.g. for the left interface (n even), for \E n \ < A we find 
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In addition to the occupation functions f n ,fn we intro- 
duce internal counting fields A„ , A„ which are Lagrange 
multipliers preventing charge accumulation inside the 
normal region. The action Q has to be varied over all 
possible configurations of /„ 's and A n 's ^3 • In the semi- 
classical regime, we may calculate the action in the saddle 
point approximation using the equations of motion |l4j , 
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These equations form an infinite system of coupled non- 
linear equations which have to be solved for arbitrary x 
[the occupation of incoming quasiparticles is unity (zero) 
at E < —A (E > A)]. The solutions are then sub- 
stituted back into Eq. to obtain S(x)- The Fourier 
transform J5J yielding the probability distribution P(Q) 
can be carried out in the stationary phase approxima- 
tion. The relation to the OBTK equations becomes clear 
if we calculate the average current / = dS(x)/dx\x=o/ T - 
For this purpose, it is sufficient to solve the saddle point 
equations for x = 0, giving A n , A„ = 0. The first two 
derivatives in Eq. © then become exactly the OBTK 
equations for the occupation functions. For the noise 
Si = d 2 S(x)/dx 2 \x=o/ T , this procedure is no longer ad- 
equate and the saddle point equations then have to be 
solved to first order in x- 

Numerical Results - In general, the FCS has to be 
calculated numerically. A typical distribution function 
P{Q) is shown in Fig. [5] the FCS is asymmetric, bounded 
from the left side, Q > (since charge may not flow 
against the bias) and exhibits a long tail on the right 
side (since the overall probability to undergo MAR and 
to climb up the energy ladder is small). Saddle point 
solutions for the /'s belonging to different points of the 
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FIG. 2: The FCS for V = 0.34A and T = 0.5. Main panel: 
The left- and rightgoing distribution functions at the right 
interface for a transmitted charge Q < (Q), Q = (Q) and 
Q > (Q)- Upper inset: The probability distribution P(Q) of 
transmitted charge. 




FIG. 3: Differential full counting statistics as function of the 
voltage for the MAR-system in the tunneling limit, T = 0.1. 



probability distribution show characteristic structures at 
E = ±A + nV due to the cusps in the Andreev re- 
flection probability (sec Fig. [TJ. The occupation func- 
tions / , / can be interpreted as the most probable 
configurations under the condition that a certain charge 
Q oc J dE( f — /) is transmitted. If the output is the 
mean charge Q = {Q}, the solution is exactly the one 
found by the OBTK equations. For a small amount of 
transmitted charge Q < (Q), the difference / — / be- 
comes tiny (no net current). In the opposite limit Q > 
(Q), this difference is large. 

The differential FCS, shown in Fig. [3J show a pro- 
nounced subharmonic gap structure, with cusps at volt- 
ages V = 2A/n. Note that the derivative of the probabil- 
ity distribution with respect to voltage is always negative, 
i.e. the distribution gets broader towards the low voltage 
limit and acquires a strong tail for Q > (Q). The first 
three differential cumulants d{Q n )/dV are plotted in Fig. 
2] for transparencies T = 0.1,0.5 and 0.9. All cumulants 
show a subharmonic gap structure, more pronounced for 
the noise and third cumulant than for the current [23| . 

Voltage limits - For high voltages V> A, only quasi- 
particles traversing the junction once contribute to the 
differential transport statistics. The saddle point equa- 
tions in JBJ) can then be solved analytically, giving the 
FCS [q = Q/(MVt/{2tt))} 

dlnP(Q) Mr f (l-r)(l + g) 4(1 - T)g 2 

-^V— = ^Y U —q q ln I-, 2 

(7) 

the result for the double barrier system in Fig. ^ m the 
normal state j2^|. Unlike the numerical result for inter- 
mediate voltages in Fig. [3J this high voltage result scales 
linearly with V and is bounded from above (0 < q < 1). 

For low voltage V -C A, the quasiparticles injected 
below the gap traverse the junction a large number of 
times N = 2A/V ^> 1 before being emitted above the 



gap. For sufficiently strong normal backscattering at the 
NS-interfaces, T N <C T/N, the motion in energy space, 
along the ladder in Fig. ^ becomes diffusive . We may 
introduce a scale m with 1 <C m <C N below which the 
motion is ballistic and described by an energy-dependent 
conductance per rung, G = MRa (E)/(2n[l - Ra{E)}). 
On larger scales, one can apply a diffusion approximation 
fn - fn = fn and f n+m ~ f„+m(df n /dn) (and similarly 
for the A's) to Eq. JSJ. Replacing the sum in Eq. by an 
integral (putting dn — dE/V), we obtain in analogy to 
the procedure for a normal diffusive wire fl5j the action 
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Changing variables A t—> A + xEjV and dE ~ Gdy maps 
Eq. |JSJ| onto the action of a normal diffusive wire hj with 
renormalized, voltage dependent charge e — > e(2A/V). 
This gives directly the generating function in Eq. Q 
with G = A[J A A dE G- 1 }- 1 = 3A/r 2 /(27r[16(l - T)}). 
Consequently, for low voltages, the current I = VG is 
proportional to voltage, while the noise Si = 4AG/3 
saturates a a constant value. All higher cumulants p > 3 
diverge as V 2 ~ p . 

Generalization - Importantly, the stochastic path in- 
tegral method presented here can be employed to any 
semiclassical, incoherent mesoscopic SNS-junction. In 
particular, Eq. Q gives the low voltage FCS of any SNS- 
junction with sufficient normal back-scattering. We em- 
phasize that the FCS is characterized by a single pa- 
rameter, the low-voltage conductance G. For a diffusive 
normal region with conductance Gw we find for instance 
G = [4tt(8(1 - T) + 3r 2 )/(3A/r 2 ) + 2/Gw}- 1 . For neg- 
ligible interface resistance, M/(2n) ^> G w . Eq. Q then 
extends the shot-noise result of Refs. to the FCS. 

The opposite limit Mj (2w) <C Gw corresponds to an 
incoherent chaotic cavity. One should however keep in 
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FIG. 4: The first three differential cumulants of the FCS, 
normalized with the normal state conductance, as a func- 
tion of voltage for three different transparencies T of the NS- 
interfaces. 



mind that eventually, for sufficiently low voltage, inelas- 
tic electron-electron and electron-phonon scattering will 
dominate the transport and the low-voltage result in Eq. 
will not apply. 

To access the FCS for the full voltage range, the nu- 
merical scheme presented for the OTBK-system has to 
be appropriately modified. E.g, for the normal region 
being a chaotic cavity, the only modification is to con- 
sider an isotropic occupation function /„ = /„ = /„ (and 
A„ = A„ = A„) in Eq. (JjjJ. For a diffusive SNS-junction, 
one can directly extend the semiclassical approach for 
current and noise in Refs. H 0, to the FCS. At a 
given voltage 2A/(n+l) < V < 2A/n, quasiparticles are 
transfered through the gap via n + 1 orn + 2 traversals 
across the junction, transporting n + 1 or n + 2 charges 
respectively. The transport is effectively through a dif- 
fusive wire of conductance Gw/{ n + 1) or Gw/( n + 2). 
The action is then just the sum of the actions of these two 
diffusive transport processes, with appropriate weights, 
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The two first cumulants 0, 0] / = GyyV and Si = 
2G^y(2A + V)/3 are smooth functions of voltage. How- 
ever, higher cumulants exhibit kinks at V = 2A/n. 

In conclusion, we have presented a theory for the full 
counting statistics of incoherent multiple Andrcev re- 
flection, based on the stochastic path integral approach. 
The charge is transfered in quantas of multiple electron 
charge, giving rise to a low-voltage divergence for cumu- 
lants of order three and higher. 
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